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We study the heat current J in a classical one-dimensional disordered chain with onsite pinning 
and with ends connected to stochastic thermal reservoirs at different temperatures. In the absence 
of anharmonicity all modes are localized and there is a gap in the spectrum. Consequently J 
decays exponentially with system size N. Using simulations we find that even a small amount of 
anharmonicity leads to a J ~ 1/N dependence, implying diffusive transport of energy. 

PACS numbers: 44.10.-K, 05.60.Cd, 05.70.Ln 



The effect of interactions between electrons or phonons 
on localization caused by disorder is a subject of both the- 
oretical P, [IS H 0, @, 0, B 0, M, ED, El and experimen- 
tal interest 13[. As is well known localization of eigen- 
functions or of normal modes strongly affects transport 
in materials containing random impurities. Localization, 
first discovered in electronic systems by Anderson 14j |. 



has its strongest effect in one dimensions where any fi- 
nite disorder makes all eigenstates localized [IH and one 
has an insulator. The presence of inelastic scattering, 
such as is caused by electron-phonon interactions, leads 
to hopping of electrons between localized states and gives 
rise to a finite conductivity. The question as to whether 
electron-electron interactions lead to a similar effect [H, Q 
has attracted much attention recently but is still not fully 
understood % @, 0, S, @] . 

In this paper we address the same question in the con- 
text of heat conduction by phonons and consider the ef- 
fect that phonon-phonon interactions have on localiza- 
tion. In particular we investigate the effect of anhar- 
monicities on the steady state transport of heat through 
a chain of oscillators with random masses. We focus on 
the case where the masses are subjected to an external 
pinning potential, in addition to nearest neighbor inter- 
actions. Pinning greatly enhances the difference between 
heat transport in a random chain with and without an- 
harmonicity and thus is a good testing ground for the 
effect of anharmonicity on localization. We also discuss 
the unpinned case and comment on results from some 
earlier studies. 

The Hamiltonian of our system has the form 
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where {xi,pi} denote the position and momenta of the 
particles and we set xq — xn+i = 0. The masses {mi} 
are chosen independently from some distribution p(m), 



e.g. one uniform in the interval (m — A, m + A). The 
chain is connected at its ends to two heat baths at tem- 
peratures Tt, and Tr respectively. The baths will be 
modeled by Ornstcin-Uhlcnbeck (Langevin white noise) 
reservoirs. The equations of motion of the chain are then 
given by: 
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Ixf 



h X[ 

v[(xi - Xl-l 



- k(2xi - xi-i - xi+i) 
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with r/i = T] L §i t x + r) R 5i :N , ji = 7(^4 + 5i, N ), and 
where the Gaussian noise terms satisfy the fluctuation 
dissipation relations (i]L(t)r]T,(t')) = 2 / yk,BTi J 5(t — t'), 
(1k(*)'/h(0) = 2jkBTn5(t — t'), ks being Boltzmann's 
constant. It is known that this system has a unique sta- 
tionary state [3]. 

Here we investigate the N dependence of the heat 
current and the temperature profile, in the nonequi- 
librium stationary state (NESS) of this system, when 
Tt, > Tr. The heat current from left to right is given 
by (Jjv) = T,i(f 1,1-1^1) / ( N ~ !) where / M _i is the force 
exerted by the (I — l)th particle on the Ith particle and 
(...) denotes an average over the NESS. It follows from 
stationarity that each term in the sum, I — 2, ...,7V is 
equal to (Jjv)- 

Note that Eq. © is invariant under the transformation 
T L . R -> sT L . R , {xi} -> {s 1 / 2 x l } and (/, v) -> {l, v)/s. 
This implies the scaling relation (Jn(sTl, sTr,1,v)) = 
s(Jn(Tl,Tr, si, sv)) . Thus the effect of changing non- 
linearity could be equivalently studied by changing tem- 
peratures. 



In the harmonic case, A = 
gives (Jjv) proportional to (Tl 



v = 0, the scaling 
- Tr). In this case 
the quadratic Hamiltonian can be written in the form 
H = (1/2)[P T M- 1 P + X T $X] using matrix notation 
for the mass and force matrices. The stationary heat cur- 
rent, for any given realization of disorder is then given 



2 



by the following expression [l7|, [IE 



(Jn) 



where 



k B (T L ~T R ) 



(3) 



T N (w) = 7 2 ^ 2 |Giw| 2 , 



nvcs 



with the matrix G = [-lj 2 M + $ — E] 1 and = 
ijuiSimlSn + Si N }. 

The NESS in the harmonic chain without random- 
ness, i.e mi = m can be solved exactly. 
(Jn) — ► cfcs(Tx — Tr), with c > 0, for TV — > oo 
i.e. the conductivity grows linearly with N. 

The nature of (Jn) for the random harmonic chain 
without pinning was analyzed in much detail in 17], LL8 



It 

19. 
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In 



17| it was proved that, for an arbitrary 
nontrivial random mass distribution, (Jn) ~ > 0, with 



(Jn) > N- z ' 2 . The results in [17J, ll& |21 



22j are based 



on the use of the Furstenberg theorem [23| first intro- 
duced into this problem by Matsuda and Ishii [24|. (In 
fact, Casher and Lebowitz use the (Jn) — » as N — > oo 
result to prove the absence of an absolutely continuous 
part of the spectrum for the random semi-infinite chain.) 
The fact that the random ID system has no extended 
states is of course well-known. The mathematical proof 
given in [25| shows that such a chain has pure point spec- 
trum, i.e. all the eigenfunctions are square integrable. 
The flux (Jn) hi this chain is carried entirely by the 
long wave-length modes with frequencies ui ~ A -1 / 2 
which "do not see" the randomness. More detailed re- 
sults about this case can be found in 17, ll 



21 



22]. In 
we be- 



particular, based on numerical evidence [1 
lieve that indeed (Jn) ~ A~ 3 / 2 . We note however that 
(Jn) m the unpinned harmonic case with disorder de- 
pends on the particular type of heat bath [ijj]. For the 
Rubin model [22j, |27| , where the baths are semi- infinite 
ordered harmonic chains in equilibrium at temperatures 
Tl(Tr) & t the left (right) ends of the system and with no 
pinning anywhere, one gets (Jn) ~ A -1 / 2 . 

Considering now the random mass harmonic pinned 
case, k a > 0, one can show that (Jn) ~ e~ cN , c > 0. 
This follows from the fact that the spectrum of this 
chain now lies entirely in the interval [(^) 1 / 2 , ( fc °+ 2fc ) 1 / 2 ] 
where m(m) are the maximum (minimum) of the masses 
{mi}. This means that there is a gap in the spec- 
trum and Furstenberg's theorem implies that Tn(u>) ~ 
g -NS(w) ^ 5(uS) > for all u> in the allowed frequency 
range. An elementary calculation shows that Tn(uj) ~ 
e~ Na , a > for w outside the spectrum (true for both 
the ordered and disordered case). Hence, in the disor- 
dered harmonic pinned case, we will have asymptotically 
for almost all realizations of disorder 



^HJn) 



c = t 



(4) 



where we can interpret £ as the largest localization length. 
A rough estimate of the value of i can be obtained using 



known results for the unpinned case at small u>. Sub- 
stituting in these results the smallest allowed frequency 
ivf n = k /m we get I = 12km/ (uj 2 ^ A 2 ) for masses cho- 
sen from a uniform distribution between [m — A, m + A] . 
This can be compared with the results from a numeri- 
cal evaluation of the integral in Eq. © for the case with 
k = k = l,m=l and A = 0.2. In Fig. (1) (A = data) 
we have plotted N[(Jn)], averaged over 100 disorder re- 
alizations, as a function of system size. From the data 
we find an exponential decay with £ w 200 while our very 
rough estimate gives i « 360. 

For the anharmonic chain there are no rigorous result 
about (Jn), even for the ordered case but the general 
expectation is that (Jn) ~ N a ~ 1 for the unpinned case. 
The actual value of a is a matter of some dispute with 
values ranging in the interval a = 1/3 - 1/2 (M 0, S3] . 
The most recent simulations found a = 1/3 [3l|. 

A number of simulations for the disordered unpinned 
case with v > 0, k Q = A = show that anharmonicity 
in general tends to destroy localization [H3 |. In a recent 
work, Li et al [32| found a transition from a = at 
small anharmonicity to a ~ 0.43 at large anharmonicity. 
They did the simulations using a deterministic model for 
the reservoirs, namely Nose-Hoover thermostats. How- 
ever this model of thermostats has been shown in j33j to 
be problematic for the harmonic case and these problems 
persist for small anharmonicity as well [34| . Preliminary 
simulations carried out for the unpinned disordered an- 
harmonic case, with white noise Langevin heat baths, do 
not find evidence of a finite v transition [35| • 

We now consider the pinned case k — k Q — 1 , A ^> 
0,^ = when the system without randomness is much 
better understood. It is generally agreed that, in the ab- 
sence of disorder, models with onsite pinning and anhar- 
monicity show regular heat transport with (Jn) ~ 1/N 
[H, [13] • This is proven rigorously for the case when a 
certain amount of stochasticity is added to the dynamics 
[3^ . This model is also closer in spirit to charge trans- 
port by hopping in random media 14 , li[ . 

In our simulations with both disorder and anhar- 
monicity we fixed the values of k, Tr, (m), r y. We then 
measure time and distance in units of ((m) /k) 1 ' 2 and 
(kBTn-f) 1 / 2 /k respectively and set k Q /k = 1, T L /T R = 2. 
The only free parameters A and A are measured in units 
of (m) and k 3 / (ksTnj) respectively (in the simulations 
we set k = k a = (m) = 7 = T R = 1). In Fig. (P) we 
plot N[< Jn >] as a function of N for a fixed disorder 
strength A = 0.2 and different values of anharmonicity 
A = 0.004 — 1.0. As can be seen from our data, there is 
a dramatic increase in the heat current on introduction 
of a small amount of anharmonicity and the system-size 
dependence goes from exponential decay to a 1/N de- 
pendence implying diffusive transport. For smaller A the 
diffusive regime sets in at larger length scales. Similar 
results are obtained for the case with A = 0, v > and 
are shown in the inset of Fig. |T]). 



The simulations were done by the velocity Verlet algo- 
rithm adapted for Langevin dynamics 39] . Equilibration 
times ranged from 10 s — 2 x 10 s time steps of step-size 
0.005 — 0.01 and steady state averages were taken over 
another 10 8 — 8 x 10 s time steps. Equilibration times in- 
crease rapidly with decreasing A and with increasing N 
and these cases required the longest runs. The error bars 
correspond to sample to sample fluctuations (other errors 
are smaller) and the number of samples varied from 10 
for the small sizes to 3 for the largest (N — 4096). We 
find that the sample-sample fluctuations become smaller 
with increasing N. A measure of the relative strengths of 
anharmonicity and disorder is obtained by looking at the 
ratio of the energy scales E a = l(x A ) /4 and Ed = TA/m. 
For our parameters we estimate e = E a /Ed ss 0.3AT/A 
and for T = (T L + T R )/2 = 1.5 this gives e w 0.008 for 
A = 0.004. 

It follows from the scaling relation described earlier 
that (J N ) = (T L - T R )f(XT L ,lT R ). Hence the conduc- 
tivity, for N — > oo and Tl — > T R = T, should depend 
only on AT, k = k(IT). For T L > T R the tempera- 
ture and hence the conductivity varies across the chain, 
we therefore measure the effective thermal conductiv- 
ity K e ff = N[< Jn >], evaluated at large N (bound- 
ary temperature jumps are negligible). For small A we 
find K e ff ~ (lT) a with a < 1. Assuming that this is 
also the behaviour of k this implies, using the fact that 
KdT/dy =constant, that the temperature profile is given 
by T(y) = [T£ +1 (l - y) + Tg +1 2/]V(«+i) where y is the 
distance, scaled by N, from the left end of the chain. We 
now look at the local temperature profiles in the NESS 
obtained in the simulations. The local temperature at 
the ith site is defined by Tj = mi{xj). The temperature 
profiles for / = 0.1 for a single sample, with and without 
disorder, are plotted in Fig. ((2]). The plots (a),(b) corre- 
spond respectively to averaging over 8 x 10 8 and 32 x 10 8 
time steps with dt = 0.005 and show that the profiles 
are reasonably converged. The noise amplitude does not 
decrease much on increased averaging. We find that the 
temperature profile, for the disordered case, is consistent 
with the predicted form with a ~ 1/2 (see however be- 
low). 

For the case without disorder we find a temperature 
profile consistent with a 1/T 2 (thus a = —2) dependence 
for the conductivity. The k(T) ~ 1/T 2 dependence has 
been predicted from recent kinetic theory calculations 
and should apply for small nonlinearity [371 ] . The inset 
in Fig. ([2|) shows temperature profiles for the disordered 
chain with a smaller and a larger value of A. We find that 
these profiles have some amount of structure, possibly re- 
flecting the local mass distribution, and vary appreciably 
from sample to sample. 

Discussion: We have studied heat conduction in a dis- 
ordered pinned anharmonic chain. We find that introduc- 
tion of a small amount of phonon-phonon interactions in 
the disordered harmonic chain leads to diffusive energy 



A 

Z 

V 

z 



0.01 




FIG. 1: Plot of the disordered averaged heat current [(Jjv)] 
multiplied by iV as a function of N for different values of A. 
The inset shows results obtained for the case with interparticle 
anharmonicity. 
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FIG. 2: The temperature profile for an anharmonic chain with 
A = 0.1 and iV = 2048 is shown for both the ordered and 
disordered cases. The inset shows the temperature profiles 
for a smaller value of anharmoncity A = 0.02 (upper curve) 
and a larger value A = 0.5. In all disordered cases A = 0.2. 



transfer, i.e. the insulating chain becomes a normal heat 
conductor. We do not find evidence of the existence of 
a finite critical value of anharmonicity required for this 
transition. For small values of anharmonicity it is nec- 
essary to go to larger system sizes to see the transition 
from insulating to diffusive. Hence a transition to a lo- 
calized phase at a very small value of anharmonicity is 
possible and would be difficult to observe in simulations. 
Even assuming that no such transition occurs, the limit- 
ing behavior of k(XT) for (AT) — > cannot be obtained 
from our data. A. Pal and D. Huse have investi- 
gated this type of question for a chain of classical spins. 
Their system corresponds very roughly to the case where 
k = A = and so the interaction between particles on 
different sites only comes via v: they find a diffusivity 
which behaves like e~ c / u . We remark that for the or- 
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dered harmonic chain letting k — * gives (Jjv) ~ k for 
the unpinned case and (Jn) ~ fc 2 for the pinned case. 

We note that in the ordered chain nonlinearity leads 
to scattering of otherwise freely moving phonons thereby 
changing the system from a super-heat conductor to a 
normal one. By contrast, in the disordered chain the 
nonlinearity appears to produce effective extended states 
which lead to diffusive transfer of energy between local 
regions. How exactly this occurs is not clear. 

The diffusive transport observed in our disordered 
open system seems to differ from recent numerical stud- 
ies of spreading of localized energy pulses in disordered 
nonlinear lattices which suggest sub-diffusion [gj or even 
absence of diffusion [9]. To understand this we have 
also performed some studies on our system without heat 
baths. From these we find that while a heat pulse intro- 
duced at one end of a long chain does indeed not prop- 
agate, a periodic driving at one end leads to significant 
energy transmission. This suggests that an isolated pulse 
behaves differently from a continuously driven or open 
system j4lj . 

Finally we note that experimental measurements of 
heat conduction in one-dimensional systems are now be- 
coming possible [H, Hi], HI] . Interesting measurable ef- 
fects which our study suggests are the effect of temper- 
ature and system-size on heat conduction in disordered 
wires. Both of these enhance the effect of nonlinearity 
and lead to diffusive transport. 
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